Geometric phases and quantum phase transitions 
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Quantum phase transition is one of the main interests in the field of condensed matter physics, while geomet- 
ric phase is a fundamental concept and has attracted considerable interest in the field of quantum mechanics. 
However, no relevant relation was recognized before recent work. In this paper, we present a review of the 
connection recently established between these two interesting fields: investigations in the geometric phase of 
the many-body systems have revealed so-called "criticality of geometric phase", in which geometric phase as- 
sociated with the many-body ground state exhibits universality, or scaling behavior in the vicinity of the critical 
point. In addition, we address the recent advances on the connection of some other geometric quantities and 
quantum phase transitions. The closed relation recently recognized between quantum phase transitions and some 
of geometric quantities may open attractive avenues and fruitful dialog between different scientific communities. 

PACS numbers: 



I. INTRODUCTION 

Quantum phase transition (QPT), which is closely asso- 
ciated with the fundamental changes that can occur in the 
macroscopic nature of matter at zero temperature due to small 
variations in a given external parameter, is certainly one of 
the major interests in condensed matter physics. Actually, the 
past decade has seen a substantial rejuvenation of interest in 
the study of quantum phase transition, driven by experiments 
on the cupric superconductors, the heavy fermion materials, 
insulator-superfluid transition in ultrocold atoms, organic con- 
ductors and related compounds 01 Q]. Quantum phase transi- 
tions are characterized by the dramatic changes in the ground 
state properties of a system driven by quantum fluctuations. 
Traditionally phases and phase transitions are described by the 
Ginzburg-Landau symmetry-breaking theory based on order 
parameters and long range correlation. Recently, substantially 
effort has been devoted to the analysis of quantum phase tran- 
sitions from other intriguing perspectives, such as topological 
orderyfl, quantum entanglementyl H] , geometric phases||5l [g] 
and some other geometric quantities!?! M l9l [Tol 

It is well-known that geometric ideas have played an im- 
portant role in physics. For example, Minkiwski's geomet- 
ric reformulation of special relativity by means of a space- 
time geometry was very useful in the construction of gen- 
eral relativity by Einstein. In this paper we will address an- 
other example: the study of quantum phase transition from 
the perspective of geometric phase (GP) factors. Actually,the 
phase factor of a wave function is the source of all interfer- 
ence phenomena and one of most fundamental concepts in 
quantum physics. The first considerable progress in this field 
is achieved by Aharonov and Bohm in 1959 lfl2ll . They pro- 
posed that the loop integral of the electromagnetic potentials 
gives an observed nonintegrable phase factor in electron inter- 
ference experiments. By using the non-Abelian phase factor, 
Yang reformulated the concept of gauge fields in an integral 
formalism in 1974 lfl3ll . and then Wu and Yang showed that the 
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gauge phase factor gives an intrinsic and complete description 
of electromagnetism. It neither underdescribes nor overde- 
scribes it lfTill . The recent considerable interests in this field 
are motivated by a pioneer work by Berry in 1984 ||T5ll , where 
he discovered that a geometric phase, in addition to the usual 
dynamical phase, is accumulated on the wave function of a 
quantum system, provided that the Hamiltonian is cyclic and 
adiabatic. It was Simon who first recognized the deep geo- 
metric meaning underlying Berry's phase. He observed that 
geometric phase is what mathematicians would call a U{1) 
holonomy in the parameter space, and the natural mathemati- 
cal context for holonomy is the theory of fiber bundles flrjll . A 
further important generalization of Berry's concept was intro- 
duced by Aharonov and Anandan lfl7ll . provided that the evo- 
lution of the state is cyclic. Besides, Samuel and Bhandari in- 
troduced a more general geometric phase in the nonadiabatic 
noncyclic evolution of the system 1 1811 . Now the applications 
of Berry phases and its generalizations II31[l6l[l7[ll8l 
can be found in many physical fields, such as optics, magnetic 
resonance, molecular and atomic physics, condensed matter 
physics and quantum computation, etc. I2lll22ll23[l24ll25[l26ll . 

Very recently, investigations in the geometric phase of the 
many-body systems have revealed so-called "criticality of ge- 
ometric phase"^ 0], in which geometric phase associated 
with the ground state exhibits universality, or scaling behav- 
ior, around the critical point^. The closed relation between 
quantum phase transitions and geometric phases may be un- 
derstood from an intuitive view: quantum phase transitions 
occur for a parameter region where the energy levels of the 
ground state and the excited state cross or have an avoided 
crossing, while geometric phase, as a measure of the curva- 
ture of Hilbert space, can reflect the energy structures and 
then can capture certain essential features of quantum phase 
transitions Ja]. 

A typical example to show the significant connection be- 
tween geometric phase and quantum phase transition is one- 
dimensional XY spin chainj5l [(J. Since the XY spin chain 
model is exactly solvable and still presents a rich structure, it 
has become a benchmark to test many new concepts. The XY 
spin chain model and the geometric phase that corresponds to 
the quantum phase transition have been analyzed in detail in 
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Ref.JH @]. The XY model is parameterized by 7 and A (see 
the definitions below Eq.©). Two distinct critical regions 
appear in parameter space: the segment (7, A) = (0, (0, 1)) 
for the XX chain and the critical line A c = 1 for the whole 
family of the XY model ifll [27ll . It has been shown that geo- 
metric phase can be used to characterize the above two crit- 
ical regions JH 0, HU. As for the first critical region, a non- 
contractible geometric phase itself f(i [28]] or its difference be- 
tween the ground state and the first excited state j5|] exists in 
the XX chain if and only if the closed evolution path circulates 
a region of criticality. There are much more physics in the sec- 
ond critical region since second order quantum phase transi- 
tion occur there. The geometric phase of the ground state has 
been shown to have scaling behavior near the critical point of 
the XY model. In particular, it has been found that the geo- 
metric phase is non-analytical and its derivative with respect 
to the field strength A diverges logarithmically near the crit- 
ical line described by A c = 1. Together with a logarithmic 
divergence of the derivative as a function of system size, the 
critical exponents are derived based on the scaling ansatz in 
the case of logarithmic divergence |29|| . Furthermore, univer- 
sality in the critical properties of geometric phase for a family 
of XY models is verified. These results show that the key in- 
gredients of quantum criticality are present in the ground-state 
geometric phase and therefore are indicators of criticality of 
geometric phase @|. 

(a) 
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FIG. 1: Schematic diagrams of the physical patterns reviewed in the 
paper, (a) Pattern I: N spins in one-dimensional chain is the whole 
system. The geometric phase of the whole N-spin system has close 
relation with quantum phase transitions of the whole system, (b) Pat- 
tern II: N spins are arranged in a circle and a test qubit in the center 
possesses homogeneous coupling with all N-spin in the ring. The 
geometric phase of the test qubit may be used to locate the criticality 
of quantum phase transition exhibits in the N-spin system. Depend- 
ing on the couplings between the spins, N-spin chain (ring) in (a) 
and (b) can be classified as the XY model, the Dicke model and the 
Lipkin-Meshkoc-Glick model. All these three models exhibit quan- 
tum phase transitions which features can be captured by the geomet- 
ric phases or some other geometric quantities. 

Motivated by these results in the XY modelJUl^], the crit- 
icality of geometric phase for other many-body models are 
investigated^ [Hi!! [H0H!- Roughly speaking, there 
are two patterns (see Fig.l) in literature to investigate the 
criticality of geometric phase in the many-body systems: (i) 



Pattern I is order to investigate the relation between geo- 
metric phase of the whole many-body system and the sys- 
tem's quantum phase transition. As illustrate in Fig.l (a), 
the geometric phase of the whole N-spin system is calculated 
and its scaling features in the vicinity of critical points are 
discussedJl U HI US HJ [H. (b) Pattern II is concerned 
with the geometric phase of a test qubit as shown in Fig. 1(b). 
N spins are arranged in a circle and a test qubit in the cen- 
ter possesses homogeneous coupling with all N-spin in the 
ring lfTll l33i l34tl . The geometric phase of the test qubit may 
be used to locate the criticality of quantum phase transition 
exhibiting in the N-spin system ic l34ll . Depending on the 
couplings between the spins, N-spin chain (ring) in (a) and (b) 
can be classified as the XY model, the Dicke model and the 
Lipkin-Meshkoc-Glick model. All these three models exhibit 
quantum phase transitions, whose features can be captured by 
the geometric phases in both patterns I and II. 

Furthermore, the study of QPTs by using other geometric 
quantities, such as quantum overlap (quantum fidelity)||7|, the 
Riemannian tensor[8|] etc., has been put forward and fruitful 
results have been reported in literature. In particular, GP is 
a imagine part of quantum geometric tensor and quantum fi- 
delity is a real part, therefore a unified theory of study QPTs 
from the perspective of quantum geometric tensor has been 
developed lfloll . 

In this paper we will review some aspects of the theoretical 
understanding that has emerged over the past several years 
towards understanding the close relation between GPs and 
QPTs. In section 2, we present the connection between Berry 
curvature and QPs. Section 3 describes the detailed relation 
between QPT and GP in the patter I. Section 4 discusses the 
results in the patter II. Finally, Section 5 presents some dis- 
cussion and perspective in the topic reviewed in this paper, in 
particular, we address the recent advances in the connection 
of some other geometric quantities and QPTs. 



II. BERRY CURVATURE AND QUANTUM PHASE 
TRANSITIONS 

Let us first address the close relation between quantum 
phase transitions and geometric phases from an intuitive 
view. Consider a generic many-body system described by the 
Hamiltonian H(rj) with r\ a dimensionless coupling constant. 
For any reasonable 77, all observable properties of the ground 
state of H will vary smoothly as 77 is varied. However, there 
may be special points denoted as i] c , where there is a non- 
analyticity in some properties of the ground state at zero tem- 
perature, 7] c is identified as the position of a quantum phase 
transition. Non-analytical behavior generally occur at level 
crossings or avoided level crossings Q|]. Surprisingly, the ge- 
ometric phase is able to capture such kinds of level structures 
and is therefore expected to signal the presence of quantum 
phase transitions. To address this relation in greater detail, 
we review geometric phases in a generic many-body system 
where the Hamiltonian can be changed by varying the param- 
eters R on which it depends. The state \ip(t)} of the system 
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evolves according to Schrodinger equation 

ihd t \m) = H(R(t))\m)- 



(i) 



At any instant, the natural basis consists of the eigenstates 
|n(R)) of H(R) forR = R(i), that satisfy F(R)|n(R)) = 
J5 n (R)|ra(R)) with energy E n (R) (n = 1,2,3- ••)■ Berry 
showed that the GP for a specific eigenstate, such as the 
ground state (\g) = |1)) of a many-body system we con- 
cern here, adiabatically undergoing a closed path in parameter 
space denoted by C, is given by [ 15] 



P a {G) = -J f V g (R) ■ dS, 



(2) 



where dS denotes area element in R space and V g (R) is the 
Berry curvature given by 



V g (R)=Im^ 



(g\V R H\n)(n\V R H\g) 



(3) 



The energy denominators in Eq.fO show that the Berry curva- 
ture usually diverges at the point of parameter space where 
energy levels are cross and may have maximum values at 
avoided level crossings. Thus level crossings or avoided level 
crossings (seem Fig. 2), the two specific level structures re- 
lated to quantum phase transitions, are reflected in the geom- 
etry of the Hilbert space of the system and can be captured by 
the Berry curvature of the ground state. However, although the 
Berry curvature is gauge invariant and is therefore an observ- 
able quantity, no feasible experimental setup has been pro- 
posed to directly observe it. On the other hand, the area in- 
tegral of Berry curvature, i.e., the geometric phase may be 
measured by the interference experiments. Therefore, rather 
than the Berry curvature, hereafter we will focus on the rela- 
tion between geometric phase and quantum phase transition, 
and therefore the proposed relation between them may be ex- 
perimentally tested. 
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FIG. 2: Schematic representation of the energy level for the many- 
body systems. The energy levels of the ground state and the excited 
state cross in (a) and have an avoided crossing in (b). On the one 
hand, quantum phase transition occurs at level crossings or avoided 
level crossings, which represents by the parameter rj c ; on the other 
hand, the Berry curvature usually diverges or may have maximum 
values at the point of parameter r] c . 



III. PATTERN I: QPT AND GP OF THE MANY-BODY 
SYSTEMS 



In this section we review the closed relation between QPTs 
and GPs for the Pattern I, as shown in Fig.l (a), where the N- 
spin chain can be classified as the XY model, the Dicke model 
and the Lipkin-Meshkoc-Glick model. 



A. The XY spin chain 

Our first example is one-dimensional XY spin chain investi- 
gated in detail in Ref.|0] . The XY model concerns N spin- 1/2 
particles (qubits) with nearest neighbor interactions and an ex- 
ternal magnetic field. The Hamiltonian of the XY spin chain 
has the following form 



H 



M 

£ 

j=-M 



1 + 7 



7 



3+1 



'3 "3 + 



(4) 

where er^ 1 (/i = x, y, x) are the Pauli matrices for the jth 
spin, 7 represents the anisotropy in the x — y plane and A 
is the intensity of the magnetic field applied in the z direction. 
We assume periodic boundary conditions for simplicity and 
choose N (= 2M + 1) odd to avoid the subtleties connected 
with the boundary terms. Nevertheless, the differences with 
other boundary conditions and the even N case are the order to 
0(1/N) and then negligible in the thermodynamic limit where 
quantum phase transitions occur H l27ll . This XY model en- 
compasses two other well-known spin models: it turns into 
transverse Ising chain for 7=1 and the XX (isotropic XY) 
chain in a transverse field for 7 = 0. 

In order to derive the geometric phase of ground state in this 
system, we introduce a new family of Hamiltonians that can 
be described by applying a rotation of cf) around the z direction 
to each spin |aJ, i.e., 



M 



H, 



4, = Uj,HU+, 



' expH^/2). (5) 



--M 



The critical behavior is independent of <fi as the spectrum 
(see below) of the system is <f> independent. This class of 
models can be diagonalized by means of the Jordan-Wigner 
transformation that maps spins to one-dimensional spinless 
fermions with creation and annihilation operators Oj and a] 

via the relations, a 3 = (IIi<j a i^ a ] d ^3l- Due t0 the 
(quasi) translational symmetry of the system we may intro- 
duce Fourier transforms of the fermionic operator described 
by d k = -^L J2j aj exp(-i2-Kjk/N) with k = —M, • • • , M. 
The Hamiltonian iJ^ can be diagonalized by transforming 
the fermion operators in momentum space and then using the 
standard Bogoliubov transformation. In this way, we obtain 
the following diagonalized form of the Hamiltonian, 



ff = j> fc (4cfc-i), 



(6) 
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where the energy of one particle excitation is given by 

A fe = \J(\- cos(2nk/N)) 2 + 7 2 sin 2 (27rfc/iV) (7) 

and Cfc = dk cos ^ — id}_ k e 21 ^ sin 4r with the angle 9k de- 
fined by cos 9k = (cos 2^ - A)/A fc . 

The ground state \g) of 77^ is the vacuum of the fermionic 
modes described by Ck\g) = 0. Substituting the operator Ck 
into this equation, one obtains the ground state as 

\g) = f] (cos^|0) fc |0)_ fc - le 2 ^sm^|l) fc |l)_ fc ) , 
fc=i ^ ' 

(8) 

where |0)fc and \l)k are the vacuum and single excitation of 
the fcth mode, respectively. The ground state is a tensor prod- 
uct of states, each lying in the two-dimensional Hilbert space 
spanned by |0)fc|0)_fc and The geometric phase of 

the ground state, accumulated by varying the angle <f> from 
to ir (Because the Hamiltonian has bilinear form, is it 
periodic in <j> ), is described by 

Pg = -jj j (</!«#. (9) 
The direct calculation shows JH 

M 

^ = ^E( 1 - cos ^)- (10 > 

k=l 

The term = 7f(l — cos 9k) is a geometric phase for the fcth 
mode, and represents the area in the parameter space (which is 
the Bloch sphere) enclosed by the loop determined by (9k , 4>)- 
To study the quantum criticality, we are interested in the ther- 
modynamic limit when the spin lattice number N — > oo. In 
this case the summation jj J2k=i can ^ e re pl ace d by the in- 
tegral i J* dip with <p = 2^ ; and then the geometric phase 
in the thermodynamic limit is given by 

P 9 = j (l-coB9 v )dip, (11) 
Jo 

where cosf?^ = (cos ip — A)/A y with the energy spectrum 

Ay, = \J(\ — cos ip) 2 + 7 2 sin 2 ip. 

As for quantum criticality in the XY model, there are two 
regions of criticality, defined by the existence of gapless ex- 
citations in the parameter space (7, A): (i) the XX region of 
criticality described by the segment (7, A) = (0, (0, 1)); (ii) 
the critical line A c = 1 for the whole family of the XY model. 
For the second critical region, we need to distinguish two uni- 
versality classes depending on the anisotropy 7. The critical 
features are characterized in term of a critical exponent v de- 
fined by £ ~ I A — A c | _l/ with £ representing the correlation 
length. For any value of 7, quantum criticality occurs at a crit- 
ical magnetic field A c = 1. For the interval < 7 < 1 the 
models belong to the Ising universality class characterized by 
the critical exponent v = 1, while for 7 = the model be- 
longs to the XX universality class with v = 1/2 JH|27t]. The 



close relation between geometric phase and quantum critical- 
ity for the first region has been addressed in Refs.Jll l28ll . 
here we mainly review the results for the second region, which 
is clearly more interesting in the sense that the second order 
quantum phase transitions occur there. 




FIG. 3: (color online), (a) Geometric phase j3 g of the ground state (b) 
and its derivative dj3 g /d\ as a function of the Hamiltonian parame- 
ters A and 7. The lattice size N — 10001. There are clear anomalies 
for the derivative of geometric phase along the critical line A c = 1. 

To demonstrate the relation between geometric phase and 
quantum phase transitions, we plot geometric phase j3 g and 
its derivative dj3 g /dX with respect to the field strength A and 
7 in Fig. 3. A significate feature is notable: the nonanalytical 
property of the geometric phase along the whole critical line 
A c = 1 in the XY spin model is clearly shown by anomalies 
for the derivative of geometric phase along the same line. 




Q.T O.S 0.9 1.0 1.1 1.2 1.3 



FIG. 4: (color online). The derivatives dj3 g /d\ for the Ising model 
(7 = 1) as a function of the Hamiltonian parameter A. The curves 
correspond to different lattice sizes N — 21, 101, 501, 1001, 00. 
With increasing the system sizes, the maximum becomes more pro- 
nounced. The inset shows that the position of the maximum changes 
and tends as TV" 1,803 towards the critical point A c = 1. 

To further understand the relation between geometric phase 
and quantum criticality, we study the scaling behavior of geo- 
metric phases by the finite size scaling approachj29|. We first 
look at the Ising model. The derivatives d(3 g /d\ for 7 = 1 
and different lattice sizes are plotted in Figj4] There is no real 
divergence for finite N, but the curves exhibit marked anoma- 
lies and the height of which increases with lattice size. The 
position X m of the peak can be regarded as a pseudo-critical 
point l29ll which changes and tends as J\f~ lm3 towards the 
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critical point and clearly approaches A c as N — > oo. In addi- 
tion, as shown in Ref [6], the value of d(3 g /d\ at the point A m 
diverges logarithmically with increasing lattice size as: 



-^k ra « KilniV + const., 
d\ 



(12) 



with K\ = 0.3121. On the other hand, the singular behavior 
of dPg/dX for the infinite Ising chain can be analyzed in the 
vicinity of the quantum criticality, and we find the asymptotic 
behavior as 



d\ 



«2 In | A — A c 



const. 



(13) 



with K2 = —0.3123. Accordin g to the scaling ansatz in the 
case of logarithmic divergence B29I1 . the ratio |«;2/ki| gives 
the exponent v that governs the divergence of the correlation 
length. Therefore, v ~ 1 is obtained in our numerical calcu- 
lation for the Ising chain, in agreement with the well-known 
solution of the Ising model Il27ll . 

A cornerstone of QPTs is a universality principle in which 
the critical behavior depends only on the dimension of the sys- 
tem and the symmetry of the order parameter. The XY model 
for the interval 7 E (0, 1] belong to the same universality class 
with critical exponent v = 1. To verify the universality prin- 
ciple in this model, the scaling behavior for different values 
of the parameter 7 has been numerically calculated in Ref. 
|0]. The results there shown that the asymptotic behaviors are 
still described by Eqs. (fT2l and ( fT3l l with K\ and k 2 being 
7-dependent constants, and the same critical exponent v = 1 
can be obtained for any 7 £ (0, 1]. 

Comparing with the 7 7^ case, the nature of the diver- 
gence of d/3 g /d\ at the critical point (7 = 0, A = 1) belongs 
to a different universality class, and the scaling behavior of 
geometric phase can be directly extracted from the explicit ex- 
pression of the geometric phase in the thermodynamic limit. 
The geometric phase under the thermodynamic limit can be 
obtained explicitly from Eq.dTTb for 7 = as 
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2tt, (A < 1) 

2ir - 2arccos(A), (A > 1) 



(14) 



However, it appears from Eq.dTOb that the geometric phase j3 g 
is always trivial for strictly 7 = and every finite lattice size 
M, since 9k = or tt for every k. The difference between 
the finite and infinite lattice sizes can be understood from the 
two limits N — > 00 and 7^0. Assume 7 = e with e an 
arbitrary small but still finite value, then we can still find a 
solution ipo (it implies N — > 00) for cos tpo — A = but A Vo = 
eyl — A 2 ^ for A ^ 1. Then a tt geometric phase appears 
for such tpo since 9 ipa = tt/2. Since df3 g /d\ = V2(l — A) -1 / 2 
(A — > 1~), we can infer the known result that the critical 
exponent v —1/2 for the XX model. 

Furthermore, we can confirm the known equivalent zv = 1 
between v and the dynamical exponent z from the calculations 
of geometric phases. The dynamical behavior is determined 
by the expansion of the energy spectrum, i.e., A^^o ~ <f z [1 + 
(ip£)~ z ]. Then z = 1 for 7 g (0, 1] and z — 2 for 7 = are 



found by the expansion of A v in the case tp — > 0. So we have 
zv = 1, which is indeed the case for the XY criticality Ql. 

Therefore, the above results clearly show that all the key 
ingredients of the quantum criticality are present in the geo- 
metric phases of the ground state in the XY spin model. 



B. The Dicke model 

Our second example is the Dicke model |36jl studied in 
Ref. lf30l [3lll . It consists of N two-level (qubit) systems cou- 
pled to a single Bosonic mode. The Hamiltonian is given by 

(ft=l) 



A 



H = ua + a + AJ X + —i={a ] + a)J : 



(15) 



where a, a + are the annihilation and creation operators of the 

Bosonic mode, respectively; J X:Z = Ylj=i a x z w ^ a i z 
ing the Pauli matrices for the qubit j are collective angular 
momentum operators for all qubits; A denotes the coupling 
strength between the atom and field; The parameters A and 
u> represent the transition frequency of the atom and Bosonic 
mode frequency, respectively. The prefactor 1 / ^/N is inserted 
to have a finite free energy per atom in the thermodynamic al 
limit N — * 00. This Hamiltonian is canonically equivalent to 
the Dicke Hamiltonian by a t t/2 rotation around the y axis. 

As illustrated in Ref s. lf37l [38tl . exact solutions may be ob- 
tained in the thermodynamic limit by employing a Holstein- 
Primakoff transformation of the angular momentum algebra. 
In the thermodynamical limit, the Dicke Hamiltonian under- 
goes a second quantum phase transition at the critical point 
A c = y/wA/2. When A < A c , the system is in its normal 
phase in which the ground state is highly unexcited, while 
A > A c , the system is in its superradiant phase in which both 
the bosonic field occupation and the spin magnetization ac- 
quire macroscopic values. 

Similarly to the XY spin model, in order to investigate the 
geometric phase one changes the original Hamiltonian by the 
unitary transformation TJ$ — exp(—i<fiJ x /2) where </> is a 
slowly varying parameter, and then the transformed Hamil- 
tonian is given by 

H+ = UlHU+ = |b 2 + 9 2 +B-J], (16) 

where the Hamiltonian of the free bosonic field is expressed 
in terms of canonical variables q = (aJ + a)/\f2 and p — 
i(a if — a)/\[2 that obey the standard quantization condition 



B 



(A 



Lq 



Lq 



-) with dimensionless 



VN sin 4> ' s/N cos </> ' 

parameters D = 2A/uj and L — 2^/2\/ uj is an effective 
magnetic field felt by the qubits. 

In the adiabatic limit, the geometric phase associated with 
the ground state of the system can be obtained by the Born- 
Oppenheimer approximation lf3QL [39h . In this case, the total 
wave function of the ground state of the system can be ap- 
proximated by 



(17) 



\^tot) = / dqip{q)\q) <g) \x(q,<t>))- 
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Here the state \x(q, <p) is the state of the adiabatic equation of 
the qubit ("fast") part for each fixed value of the slow variable 
q, i.e., 



B.J\x(q,<p))=E(q)\ X (q,<p)) 



(18) 



with E(q) the eigenenergy. It can be proven that the state 
\x(q, ¥>)) can be expressed as a direct product of N qubits as 
\x(q, f)) = ®f—i\x{qi l v 9 ))j: an d the state of each qubit can 
be written as 

OL OL _■ 

\x(q, V>))j = sin -| t)j - cos -e rq \ Dj 



with cos a = Lq cos (f)/(\/NE(q)) and tan 77 = 
Lq sin (f) / (\/N D) . On the other hand, the ground state 
wave function for the oscillator ip(q) is governed by one- 
dimensional time-independent Schrodinger equation 



Had\<p{q)) = 



2 \dq- 



+ q A -NE(q) =eo\<p{q)), 



where £0 is the lowerest eigenvalues of the adiabatic Hamilto- 
nian H ad . 

Once the total wave function of the ground state is derived, 
the geometric phase /3 g of the ground state may be derived by 
the standard method as j3 g = i § (iptot\d/ d(f>\iptot) dip, and the 
final result is given by 



p g = Nir[l + 



N 



In the thermodynamic limit, one can show that 



N 



N^oo 



0, 



(a<l) 
(a>l). 



(19) 



(20) 
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FIG. 5: The geometric phase 7 (= /3 g ) of the ground state and its 
derivative (inset) for the Dicke model with respect to the parameter 
a for different values of qubit iV and the parameter D — 10. The 
geometric phase increases with a, and there is a cusplike behavior in 
the thermodynamic limit at the critical transition point a = 1. 

The scaled geometric phase f3 g /N and its derivative with 
respect to the parameter a for D = 10 is shown in Fig|5][ 30]. 
It is evident that the geometric phase increases with increasing 
the coupling constant at the finite qubit number N, while in 
the thermodynamic limit the geometric phase vanishes when 
a < a c and has a cusplike behavior at the critical point a — 



a c . In addition, the derivative is discontinuous at the critical 
point. These results are consistent with the expected behavior 
of the geometric phase across the critical point, and therefore 
we add another unusual example to the close relation between 
geometric phase and quantum phase transition. 



C. The Lipkin-Meshkov-Glick model 

Our third example is the Lipkin-Meshkov-Glick (LMG) 
model discussed in Ref l32ll . The LMG was first introduced 
in nuclear phvsics ll40ll . The LMG model describes a set of N 
qubits coupled to all others with a strength independent of the 
position and the nature of the elements and a magnetic field h 
in the z direction, i.e., the Hamiltonian is given by 



1 



H = -^(S 2 x + 1 S 2 y )-hS z , 



(21) 



where 7 is the anisotropy parameter. S a = Y^iLi a a/'^( a 
x,y, z) and the a a is the Pauli operator, N is the total parti- 
cle number in this system. The prefactor 1/N is essential to 
ensure the convergence of the free energy per spin in the ther- 
modynamic limit. As widely discussed in the literature ( see, 
e.g., Ref. ll4lll ). this system displays a second-order quantum 
phase transition at the critical point h = 1. 

The diagonalization of the LMG Hamiltonian and deriva- 
tion of the geometric phase can be obtained by a standard pro- 
cedure, which can be summarized in the following steps fl32ll : 
(i) perform a rotation of the spin operators around the y di- 
rection, that makes the z axis along the so-called semiclassi- 
cal magnetization l42tl in which the Hamiltonian described in 
Eqf2T]has the minimal value in the semiclassical approxima- 
tion, (ii) Similar to the XY model and the Dicke model, to 
introduce a geometric phase of the ground state, we consider 
a system which has a rotation U(<fi) = e^ l< ^ Sz around the 
new z direction, and then the Hamiltonian becomes H((j>) = 
U> (4>)HW (4>). (iii) then we use the Holstein-Primakoff rep- 
resentation, 



SM 
§ + {<t>) 



N/2 - a) a, 
(N 



a) e 



(22) 



in which a) is bosonic operator. Since the z axis is along 
the semiclassical magnetization, a'a/N 1 is a reason- 
able assumption under low-energy approximation, in which 
N is large but finite, (iv) the Bogoliubov transformation, 
which defines the bosonic operator as b(<j>) — coshxae 1 ^ + 

2r/A with A = sin 2 6 - 

2 n 

— . These procedures diag- 



sinhxa^e ', where tanhlx 
7 +co S 2 e +fecos6>andr = 2_ 



onalize the Hamiltonian to a form 



H diag (<f>) = Nd + £ + A D 6 t (</))6(0), 



(23) 



where d 



Avl 



-i(sin 2 6»+2/icos 
and e 



tanh2a; = 2T/A. The ground state 
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FIG. 6: The geometric phase f3 g of the ground state for the LMG 
model as a function of the parameter (7, h) for N = 200. The 
divergence of (3 g is evident at the critical line h c = 1 



\g(4>)} is determined by the relation b(4>)\g((j>)) = 0. Sub- 
stituting b((f>) into the equation above, one finds the ground 
state, 



whole system may have the form 

H = Ht+H s + Hj, 



(26) 



where H t = /iB • a stands for the Hamiltonian of the test qubit 
in a general form, Hs represents the Hamiltonian of a many- 
body system which we are going to study, and Hj denotes 
the coupling between them. We assume that the quantum sys- 
tem described by Hs undergoes a quantum phase transition at 
certain critical points. It is expected that the geometric phase 
of the test qubit can be used to identify the quantum phase 
transition of the many-body system. A relatively general for- 
malism to show the close relation between geometric phase of 
the test qubit and quantum phase transition of the many body 
system has been developed in Ref . lf33h . For solidness, here 
we address a detailed example studied in Ref. l34ll . where the 
many -body system with the quantum phase transition is a XY 
spin chain, i.e., 



[N/2] 



n— 

■(-V2e 



(2n-l)!! e-^sinh: 



2n\\ e^coshx 
sinh x)\2n) , 



(24) 



where nil — n(n — 2)(n — 4) ■ ■ • and nl\ = 1 for n < 0. |n) 
is the Fock state of bosonic operator and the normalized 

constant is C 2 = V [7V/ n 21 2 sinh 2 tanh 2 ^" 1 ) x. 

The geometric phase (3 g of the ground state accumulated by 
changing <j> from to it can be derived by the standard method 
as shown before, and the final result is give by ||32ll 



a 



EL^^tanh^-^x 



(25) 



To have some basic ideas about the relation between the 
geometric phase and phase transition in the LMG model, the 
geometric phases f3 g as a function of the parameters (7, h) 
have been plotted in Fig.6 ll32ll . It is notable that the geometric 
phase (3, independent of the anisotropy, is divergent in the line 
h = 1, where the LMG model has been proven to exhibit a 
second-order phase transition ll4lll . The divergence of geomet- 
ric phase itself, rather then the derivative of geometric phase, 
shows distinguished character from the XY and Dicke mod- 
els. This difference stems from that the collective interaction 
in the LMG model, which is absent in the XY model p"^ . 

The scaling behavior of (3 g has also been studied in Ref 



0211 . A relatively simply relation (3 g w — N is obtained there. 
Furthermore the scaling is independent of 7, which means 
that for different 7, the phase transitions belong to the same 
university class. This phenomenon is different from the XY 
model, in which the isotropic and anisotropic interactions re- 
spectively belong to different university classes |@]. 



IV. PATTERN II: GP OF THE TEST QUBIT AND QPT 

In this section, we consider a test qubit coupled to a quan- 
tum many -body svstem lfTll l33l l34ll . The Hamiltonian of the 



H t = fxa z /2 + va x /2, 



(27) 



M 



1 



l=-M 



1 V V 

—°l a l+l 



H, 



N 

-y 

1=1 



a a, 



(28) 



(29) 



where the Pauli matrices a x ' v ' z and a^ v ' z denote the test 
qubit and the XY spin chain subsystems, respectively. The 
parameter 77 represents the coupling strength between the test 
qubit and all spins (qubits) in the spin chain. This model 
is similar to the Hepp-Coleman model i43ll . which was ini- 
tially proposed as a model for quantum measurement, and its 
generalizationlH El] . 

Following Ref.[15], we assume that the test qubit is ini- 
tially in a superposition state \(f> t (Q)) = c g \g) + c e |e), where 

\g) = (sin^,-cos^) T and |e) = (cos^,sin^) T with 
6 a = tan _1 (^//i) are ground and excited states of H t , re- 
spectively. The coefficients c g and c e satisfy the normaliza- 
tion condition, \c g \ 2 + |c e | 2 = 1. Then the evolution of the 
XY spin chain initially prepared in \ip(Q)), will split into two 
branches \ip a (t)) = cxp(~iH a t)\(p(0)) (a = g,e), and the 
total wave function is obtained as \ip{t)) = c g \g) (g> \ t p g (t)) + 
c e |e) (g> \<p e (t)). Here, the evolutions of the two branch wave 
functions \ip a (t)) are driven, respectively, by the two effective 
Hamiltonians 



N 



H g = {g\H\g) = H s -5j2°i-^ 



H„ 



(e\H\e) 



1=1 



N 



+ A, 



(30) 



(31) 



1=1 
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where A = y/ fx 2 + u 2 /2 and 5 = r/cosOo/N. Both H g and 
H e describe the XY model in a transverse field, but with a 
tiny difference in the field strength. Similar to the method to 
diagonalize the standard XY spin chain addressed in the patter 
I, the ground states of the Hamiltonians H a are given by 



I Go 



M 

n 

k=l 



■|0) fc |0)_ fc + zsin^|l) fc |l)_ fe 



(32) 



where cos6>£ = e£/A£ with A£ = 

£fc,Q = A - cos^ + k q S (K g = ~K e = 1). |0) fc and |l) fc 
are the vacuum and single excitation of the fcth mode, respec- 
tively. Here dj. is similarly defined as the standard XY model 
(see section 3.1). 

Now we turn to study the behaviors of the geometric phase 
for the test qubit when the XY spin chain is at its ground state. 
Due to the coupling, it is expected that the geometric phase for 
the test qubit will be profoundly influenced by the occurrence 
of quantum phase transition in spin-chain environment. Since 
we are interesting to the quantum phase transition, which is 
the property of the ground state, we assume that the XY spin 
chain is adiabatically in the ground state \G g ({9l})) of H g . 
In this case the effective mean-field Hamiltonian for the test 
qubit is given by 



H, 



eff 



= Hi 



{G^H^Gg 

M 



k=l 




(33) 
(34) 



In order to generate a geometric phase for the test qubit, as 
usual, we change the Hamiltonian by means of a unitary trans- 
formation: U ((f>) = cxp i^o^j , where is a slowly vary- 
ing parameter, changing from to it. The transformed Hamil- 
tonian can be written as H((f>) = U + (4>)H e ffU(4>), i.e., 

H{4>) = (f + ^f>s0^ ^+|(^cos0-^sin0). 

(35) 

Then the eigen-energies of the effective Hamiltonian for the 
test qubit are given by 



E, 



± 



\ 



k=l j 



v- 

T 



(36) 



and the corresponding eigenstates are given by 



\g) = 



where sin# = v/2E e . 

The accumulated ground-state geometric phase (3 g for the 
test qubit by varying (j> from zero to ir can be derived from the 
standard integral J™ {G g \d<p\G g )d(p, and it is easy to find that 



09 



1 



l i + Ai 1 f{\ 1 ,N) 
v /[m + 4 ?? /(A,7V)] 2 + ^ 



(38) 




FIG. 7: (a) Ground-state geometric phase f3 g of the test qubit and 
(b) its derivative df3 g /d\ as a function of the spin-chain parameter 
(A, 7). The anomalies for the derivative of geometric phase is clear 
along the critical line A c = 1. The other parameters: /i = 0.1, 
v = 2, and r\ = 0.5. 




FIG. 8: The derivatives dj3 g /d\ for the test qubit which is coupling 
to the Ising spin chain (7 = 1), with respect to the parameter A 
for different lattice sizes N — 12, 51, 251, 501, 00. With increasing 
the system sizes, the maximum becomes more pronounced and the 
position of the maximum clearly approaches A c = 1 as N — > 00. 
The inset shows the size scaling of the position of the peak occurred 
in dp g /dX (circles) and the function /(A, 7, N) (squares). 



where /(A, 7, N) = -h J2kLi cos ^Jf'- ^ n tne thermodynamic 
limit N — > 00, the summation in /(A, 7, N) can be replaced 
by the integral as follows: 



f(\, 1 ,N)\ N ^ 00 = ±- 
Zir 



X — cos if 



\J(\- cos^) 2 + 7 



9 • 2 
z sin (p 



-dip. 



(39) 



The geometric phase j3 g and its derivative df3 g /dX with re- 
spect to the parameter (A, 7) of the XY model are plotted in 
Fig.7. As expected, the nonanalytic behavior of the geomet- 
ric phase and the corresponding anomalies in its derivative 
(37) dP g /dX along the critical lines A c = 1 are clear. All these 



features are very similar to those in the XY spin chain in pat- 
ter I (see section 3.1). 

To further understand the relation between GPs and QPTs 
in this system, let us consider the case of XX spin model 
(7 = 0) in which geometric phase can be analytically derived. 
In the thermodynamic limit, the function /(A, 7, N) in Eql39l 
can be derived explicitly for 7 = as / = 1/2 — arccos(A) /tt 
when A < 1 and / = 1/2 when A > 1. In this case, the 
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geometric phase of the test qubit is given by 



/i-f 2g[l — 2 arccos( A) / 7r] 



V(M+2g[l-2£ 
V(m+2s) 2 +^ 5 



s(A)/^]) 2 +^ 



(A<1) 



(A>1) 
(40) 

which clearly shows a discontinuity at A = A c = 1. The 
derivative d(3 g /dX as a function of A for 7=1 and differ- 



ent lattice sizes are plotted in Fig. 8 13411 . It is notable that 
the derivative of geometric phase is peaked around the criti- 
cal point A c = 1. The amplitude of the peak is prominently 
enhanced by increasing the lattice size of the spin chain. The 
size dependent of the peak position A m for d(3 g /d\ is shown 
in the inset of Fig. 8. For comparison, the size dependence of 
the peak position in A space for the derivative df / dX are also 
shown in the inset (squires). The scaling behavior of d(3 g /dX 
and df /dX are evident in the figure. All these features are 
similar to these exhibit in the XY spin chain of the patter I. 
Therefore, we can see that QPTs of the XY spin chain are 
faithfully reflected by the behaviors of the ground-state GP 
and its derivative of the coupled test qubit. 



V. SUMMARY AND CONCLUDING REMARKS 

Quantum phase transition plays a key role in condensed 
matter physics, while the concept of geometric phase is fun- 
damental in quantum mechanics. However, no relevant re- 
lation was recognized before recent work. In this paper, we 
present a review of the connection recently established be- 
tween these two interesting fields. Phases and phase tran- 
sitions are traditionally described by the Ginzburg-Landau 
symmetry-breaking theory based on order parameters and 
long rang correlation. Recent develops offer other perspec- 
tives to understand quantum phase transitions, such as topo- 
logical order, quantum entanglement, geometric phases and 
other geometric quantities. Before conclusion, we would like 
to briefly address that, rather than geometric phase reviewed 
in this paper, the deep relationship between some other geo- 
metric quantities and quantum phase transitions has also been 
revealed. 

Quantum fidelity. Recently an approach to quantum phase 
transitions based on the concept of quantum fidelity has been 
put forward [3]- In this approach, quantum phase tran- 
sitions are characterized by investigating the properties of 
the overlap between two ground states corresponding to two 
slightly different set of parameters. The overlap between two 
states can be considered as a Hilbert-space distance, and is 
also called quantum fidelity from the perspective of quan- 
tum information. A drop of the fidelity with scaling behav- 
ior is observed in the vicinity of quantum phase transition and 
then quantitative information about critical exponents can be 
extractedJUl [46J. The physical intuition behind this relation 
is straightforward. Quantum phase transitions mark the sep- 
aration between regions of the parameter space which corre- 
spond to ground state having deeply different structural prop- 
erties. Since the fidelity is a measure of the state-state dis- 



tance, the dramatic change of the structure of the ground state 
around the quantum critical point should result in a large dis- 
tance between two ground states. The study of QPTs based 
on quantum fidelity (overlap) has been reported for several 
statistical models O, [3, |47l H 53]. In addition, the dynamic 
analogy of quantum overlap is the Loschmidt echo; it has been 
shown that the Loschmidt echo also exhibits scaling behavior 
in the vicinity of the critical point ! Ill l50u5lll . 

The Riemannian tensor. It has been shown that the fidelity 
approach can be better understood in terms of a Riemannian 
metric tensor g defined over the parameter manifold [ 8 ] . In this 
approach, the manifold of coupling constants parameterizing 
the system's Hamiltonian can be equipped with a (pseudo) 
Riemannian tensor g whose singularities correspond to the 
critical regions. 

We have presented that one can study quantum phase tran- 
sitions from the perspective of some geometric objects, such 
as geometric phase, quantum fidelity and the Riemannian ten- 
sor. Surprisingly, All these approaches share the same ori- 
gin and can be therefore unified by the concept of quan- 
tum geometric tensors. We now briefly recall the formal 
setting developed in Ref.|10j]. For each element 7/ of the 
parameter manifold M. there is an associated Hamiltonian 

H(V) = En^ E n (v)\^n(v))(Mv)\ (E n+1 > E n ), act- 
ing over a finite-dimensional state space 7i. If \^(f])) repre- 
sents the unique ground state of H(rj), then one has the map- 
ping '■ -M — > H : Tj — ► 1^(77)). In this case, one can 
define a quantum geometric tensor which is a complex her- 
mitean tensor in the parameter manifold M. given by f52ll 

Qv* = (fytfol&^o) - (5 Al *o|*o)(*o|a,*o), (41) 

where the indices fx and v denote the coordinates of M.. The 
real part of the quantum geometric tensor Q is the Rieman- 
nian metric, while the imaginary part is the curvature form 
giving rise to a geometric phase lflQIl . Similar to the heuristic 
argument that we have addressed for the singularity of Berry 
curvature in the vicinity of quantum phase transition, it has 
been shown that the quantu m g eometric tensor also obeys crit- 
ical scaling behaviorlla, M, 1 1011 . Therefore, viewing quantum 
phase transitions from the perspectives of geometric phase and 
quantum fidelity can be unified by the concept of quantum ge- 
ometric tensor. 

In conclusion, we presented a review of critic ality of geo- 
metric phase established recently, in which geometric phase 
associated with the many-body ground state exhibits univer- 
sality, or scaling behavior in the vicinity of the critical point. 
In addition, we addressed that one can investigate quantum 
phase transition from the views of some typical geometric 
quantities. The closed relation recently recognized between 
quantum phase transitions and quantum geometric tensor may 
open attractive avenues and fruitful dialog between different 
scientific communities.. 
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